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manipulation ; whilst the extent of my obligation to Professor 
Turner is very far indeed from being covered by the definite 
acknowledgments I have already made. 


Some Suggestions for the Explicit Use of Direction Cosines or 

Rectangular Coordinates in Astronomical Computations. By 

H. H. Turner, M.A., F.R.S., Savilian Professor. 

I. Until the introduction of the photographic method, the 
simple computations of astronomy could be fairly regarded as 
settled in form. The form depended essentially on the large use 
of meridian observations, which determine R.A. and N.P.D. 
independently, and thus all computations and catalogues are 
arranged in terms of R.A. and N.P.D. Tables of the planets, 
which are most conveniently expressed in latitude and longitude, 
and are so expressed in the first instance, are converted into 
terms of R.A. and N.P.D. for the use of observers. 

2. But even independently of photography some incon¬ 
veniences attending the use of these particular coordinates have 
made themselves felt. The neighbourhood of the pole is always 
a difficulty, and calls for some modification of the ordinary pro¬ 
cesses. Fabritius proposed in Ast. Each. Nos. 2072-3 to use the 
coordinates 

£c=sinp cosa y=sinp sina 

in the neighbourhood of the pole, where a, p are the R.A. and 
N.P.D. of a star ; and he showed that the use of these rectangu¬ 
lar coordinates considerably simplifies the calculation of preces¬ 
sion and nutation for polar stars. His x and y may be thus 
regarded. Draw a tangent plane to the sphere at the pole, and 
orthogonally project stars from the sphere on to this plane by 
lines parallel to the Earth’s axis ; then (x, y) are the Cartesian 
coordinates of the projection referred to axes in this plane. 

3. In connection with the measurement and reduction of 
photographic plates I have urged the advantages of using 
u standard coordinates ” which are also the Cartesian coordinates 
©f a star supposed projected on a tangent plane to the celestial 
sphere ; but in this case the projection is not orthogonal * it is 
by lines radiating from the centre of the sphere. For a plate 
with centre at the pole the u standard coordinates” would be 

£=tan p cosa rj—t&np sina 

as compared with Fabritius’s 

a?=sinp cosa ?/=sinp sina 

4. Other projections on the tangent plane might be made. 
Thus the well-known stereographic projection would be repre¬ 
sented by 

X=2 tan P cosa Y=2 tan ^ sinct 
2 2 
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and it is possible that these coordinates or yet others may be 
found useful in some department of astronomical work. But 
without following generalities further in this direction, we may 
notice that both Fabritius’s and standard coordinates are closely 
related to direction cosines. If the three direction cosines of a 
star be 

Z=sinp cosa ra=sinp sina n—cosp 5 

then we have for Fabritius coordinates 

x=l y=m 

and for standard coordinates 




l 

n 


U= 


m 

n 


Now this suggests that the more explicit use of direction cosines, 
in our computations may be attended with advantages - and the 
present paper is written with the object of calling attention to 
this point by two particular examples. x 

5. The first example is that in which we are dealing with 



positions of objects on the surface of the Sun or Moon. Let G 
be the centre of the Moon ; E, the spectator on the Earth. He 
sees any point P on the surface of the Moon projected on the 
plane AB, as at p. Now E is at a great distance, and hence EP 
is nearly parallel, though not quite, to EG, By a slight correc¬ 
tion to Gp we can find the orthogonally projected coordinates of 
P, and are thus prepared for the use of direction cosines in the 
reductions. The correction required is 

jM=P?n tan niPp 

=*yCB 2 -Cm 2 . Gp . /CE 


=°i(i 

CE l 


“(cb) i°P 


, Here GB/CE is the angular radius of the Moon in circular 
measure—about -005. The remainder of the expression hds a 
maximum value when Cm 2 ==-^CB 2 , nearly; and thus the maxi¬ 
mum value of the correction is ^ 2 

*■ \ X *005 x i 6'=2 // *5, say. 
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Hence a small table suffices to give the fraction by which we 
must diminish Op. 

If, then, we measure on the Moon, or on a photograph of the 
Moon, rectangular coordinates x' and y' from the centre of the 
disc, and expressed in terms of the Moon’s radius, and with 
s/x'z+y ' 2 argument, take out from the table the appropriate 
reducing factor, we can get x, y, two of the direction cosines of 
the point, and can then find z~\f j—x^—y 2 . 

6. These (x, y, z) are coordinates of the point with axis of $ 
directed to the centre of the Earth and axes of x and y at right 
angles to it, with some orientation depending upon the way in 
which the measures were made. Let now (X, Y, Z) be direction 
cosines of the point referred to chosen axes fixed in the Moon. 
Then by the properties of direction cosines we know that 

I=za x x 4 -a^y 4 -a^z 

y=b l x+b 2 y+b. s z 

i=c l x+c 2 y+c 3 z 

with the reciprocal relations 

x=a l £ + b l r] + c l Z 

y—a^ + btfi + CzC, 

£=<z 3 £4-&3>/ + & 3 ; 

where {a { , a 2 , ci 3 . . . c 3 ) are a system of direction cosines express¬ 
ing the relationship of the two systems. We may calculate these 
coefficients from our knowledge of the Moon’s librations y or, if 
we know the coordinates (X, Y, Z) of certain points on the 
Moon, of which we determine the (ay y t £) by measures, we can 
find the coefficients (a x . . . c 3 ) by substituting for these known 
points in the equations above given. 

7. Again, our measures may only give us x and y, affected 
with certain errors—those due to refraction, or defective cen¬ 
tring of the plate, or defective scale value, or orientation. We 
shall thus get, not x and y , but two linear functions of them 
and z, viz. 

Xq = Ax +By+ Cz 4-D 
y 0 =Ex+Fy + Gz+H 


where the coefficients AB . . . H are not altogether known, though 
certain parts of them, such as that due to refraction, may be 
calculated. Without troubling to make this calculation, how¬ 
ever, we see from the linearity of all the above relations that we 
may write 

x^aH + brj+cZ+d 
yo— e % + 9% +^ 
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so that the measures x 0 and y Q (corrected to orthogonal projection 
of course) are expressible as linear functions of the coordinates 
(£, ??, £). If, then, we know the (£, 77, £) of a certain number of 
measured points we can find these eight constants of the plate 
(#, b...h)i and then using them for any other measured point we 
can find the (£, 7?, £) from the two equations of the above form 
combined with (£, 17, 4). 

8. Of this method Mr. Saunder (see p. 184) has made trial in 
measuring some beautiful lunar photographs, kindly lent by the 
Director of the Paris Observatory (M. Loewy). He found the 
method unsatisfactory, however, because there are, in fact, at 
present not enough points on the Moon whose positions are well 
determined. Mr. Saunder has therefore devised (see p. 185) an 
elegant method of his own, which bids fair to give us excellent posi¬ 
tions of a number of points ; and then it is possible that he may 
return to this method for finding the coordinates of other points. 

9. But there is a point of detail in which the method requites 
improvement, viz. it is rather troublesome to solve the three 
equations 

a^+brf + c^=x Q — d=x x , say 
e$+fo+gZ=yo—h=y l} say 
£ 2 -H 2 + 4 2 =i 

(of which the third is a quadratic) for every point. The constants 
d and h are quickly subtracted, of course, and may be dis¬ 
regarded. To avoid the solution of the quadratic I suggest the 
following procedure :— 

But for refraction and accidental errors the constants a, b...g 
would fulfil the relations 

c& 2 + 6 2 4-c 2 = 1, e 2 4-/ 2 -j-£ 2 =i, ae + bf+cg=z o 
and then putting 

1—- a 2 — e 2 =p 2 , 1 — b~ —f 2 =q 2 , 1 —c 2 —<? 2 =r 2 , 1 
we should have 

£=ax l -\ r ey l +pz l 

rt=.bx l +fy l +q% l 

£z=Lcx l +gy l +rz l 

All these calculations could be quickly made, especially with an 
arithmometer and a table of squares. 

Owing to refraction and accidental errors these relations are 
only approximately fulfilled. But if instead of x x and y x we write 

%i+)Wi + yyi &ndy l +fix l + ry 1 . 
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where X, g, r are small, we can determine (A, //, r) so as to fulfil the 
corresponding relations exactly. We have, in fact, 

(a -f aX -f eff -f (b bX + (c + +^) 2 = 1 

.=i 

(a + a\+ejji)(e + a}x + ev)+ . . . =0 


Put 1—(a 2 + 5 2 +c 2 )=A, 1—(e 2 +/ 2 +^ 2 )=C, ae + bf+cg= r—B, 

where A, B, C are at least of the same order as A, //,, v ; and neglect 
in the first instance squares of small quantities. The equations 
become 

2\=A, 2r = C, 2/l=B, 

and the test whether this approximation is sufficient is that the 
corrected values of the coefficients should sensibly fulfil the 
above relations. If not, the process may be repeated, and new 
values of A, //, v found. But since refraction is small, one approxi¬ 
mation should suffice. Putting, then, 

a=a-\-aX-\- efjL fi=b + b\+ffA &c. 
and x< l =x l +\x l +}Ay u 2/ 2 =2/i +t*M 1 A-vy^, z 2 2 =i —xf—yf, 

so that a, /? and also x 2 , y 2 , z 2 can be found without much trouble, 
we have 

£=a%+ey 2 + 7r£ 2 

r)=fix 2 + t;y 2 + KZ 2 

%=yx 2 -Any 2 4-^2 

and thus £, rj, £ can be found simply. 

10. A second example of the explicit use of linear formulae 
of this kind is afforded by the computations of the effect of 
precession and nutation. The advantages have been partially 
pointed out by Pabritius and others. In Ast. Each., No. 3610, * 
recently published, Dr. W. Ebert considers the special applica¬ 
tion of Fabritius’s formulae to stars within 20' of the pole. But 
I venture to think that the advantages of the general applica¬ 
tion of such formulae have been hitherto overlooked. Precession 
and nutation simply change the three axes of reference, so that 
if (x, y , z) are the direction cosines of a star for one date and 
(£, rj, 4 ) those for another, we have 

Z=a x x + a 2 y-\-a z z 
n — b x x-\-b 2 yAb z z 
X—c x x-\-c^y^c z z 

where the constants a l} a 2 , a z . . . c 3 are connected by the six 
well-known relations. Since the axes of reference change slowly 

p 
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a H & 2 , c 3 are nearly unity, and the others nearly zero, for any 
moderate interval. We may write, in fact, 1 

£—aj=aoj 4- a 2 y 4- a 3 z 

v—y=h x +fy+h z 

Z—z^CyX+c^+yz 


and then the coefficients on the right are all small, and only 
approximate values of x, y, z are needed on the right. 

11, Some of the coefficients are of the first order and some of 
the second. To show their general character let us neglect the 
slow motion of the ecliptic and consider the pole of the equator 
(axis of z or £) revolving round the fixed pole of the ecliptic with 
angular velocity q. Let x, y, z he the coordinates of a star at 
epoch £== o, and £, 77, '( the coordinates at time 

The pale of the ecliptic has the same coordinates in both 
systems, viz. 

x=o y = — sin u) z=cosw 

£=o rj — — s i n <*> $4 'C— COS w , 


Thus 0 =— a 2 sin w -f a z cos w 

? —sin 6>= — & 2 sin w + 5 3 cos to 

cos o>= —c 2 sin w+c 3 cos w. 

Therefore tan u)=: a J=z J* 3 — Cs ~~i . 

®'2 b 2 ’— I C 2 


Multiply the numerators and denominators by a z , 5 3 , c 3 respec¬ 


tively and add. Each ratio is found equal to- =4 and this 

-o 3 

shows that 6 3 =c 2 . The same result is obtained if we 
multiply by a 2 , b 2 , c 2 . But from multiplying by a u by, c\ we get 

each ratio equal to 1 ; whence c l =b l tan w. It will be found 

°\ ' ' # # ; ‘ 1 ' 5 

that we can now express all the quantities in terms of c 3 . Thus 


c 2 =(c 3 — 1) cot (O 



by = C, COt 61 
63 (c 3 I^ COt 61 

6 2 2 =1-6,2-42 


and so on. To find c 3 we must refer to the spherical triangle 
formed by E, the pole of the ecliptic, and P, P 2 , the two poles of 
the equator at the different epochs. In this triangle 


<»3=cos P 1 P 2 =cos P t E cos P 2 E + sin P,E sin P 2 E cos E 
=cos 2 6>4-sin2 6> cos qt. ■ 
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Hence c 2 =sin w cos w (cos qt— 1) 

C! = sin usings 
^=003 to singtf. 

We may express the values by the following scheme :— 


£ cos qt 

r) cos 00 sin qt 

£ sin 03 sin qt 


y 

— cos o) sin qt 
cos 2 03 cos qt + sin 2 <a 
cos 03 sin 03 (cos qt— 1 ) 


— sin 03 sin qt 
cos 03 sin a> (cos qt — 1 ) 
sin 2 03 cos qt + cos 2 00. 


12. Now qt is a small quantity. In a century ^=5000" 
approximately, or, say, *025 in circular measure. Thus we may 
exhibit the approximate values of powers of qt for different 
periods as follows :— 


qt 

10 yrs. 

*0025 

5° yrs. 
*0125 

100 yrs. 

•025 

30c yrs. 

'075 

(SO* 

•OOOOO63 

•OOOI563 

•000625 

•005625 

Xgty 

•OOOOOOO 

*0000020 

•0000156 

•0004188 

(.sty 

... 

*0000000 

•0000004 

*0000316 

(sty 

... 

... 

•0000000 

*0000024 


Hence, if we are applying precessions for ten years, we can 
neglect the third power of qt; for fifty years the fourth 
power; and, indeed, up to one hundred years we can neglect the 
fourth power, which, as we shall see, is generally multiplied by a 
small fraction, 

13. For simplicity let us first neglect (qt) 3 — i.e. consider pre¬ 
cessions for moderate periods comparable with ten or twenty 
years. Then putting qt cos <*>=s, tan ; since cos <.»= # 9i > 
s is less than qt, and its powers converge more rapidly. Also 
powers of >(=0*4) converge quickly. The scheme of transfor¬ 
mation now becomes 


x 

£ I— |(i-{-r 2 )s 2 

V + s 


V z 

—s —rs 

1 — —\rs 2 

~5 rs2 i-|rV 


Since s is proportional to the first power of the time, and s 2 
to the second, we see that this scheme gives in a compendious 
form both the precessions and secular variations of x 1 y % z —viz. 
the precessions are 

in x, —( y+rz)s 
in y, -f -x.s 
in z % +xr . s 

P 2 
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and the secular variations are 

in x —\x(i +r 2 )s 2 
in y — + rz)s 2 
in % —\r(y+rz)s 2 


14. Similarly, the additional terms of the third and fourth 
orders are readily written down as follows, omitting a common 
factor, ^^(i+r 2 ), for convenience in writing :— 


x 

£ s(i+r 2 ) 

V —4 

£ —4 r 


y 

4 

-f- sr 


z 

4r 

-j-sr 
-j-s . r 2 


The terms of the third order are all less than i(qt) 3 , and only 
begin to affect the seventh place of decimals after about thirty 
years. 

The largest term of the fourth order is (qt)\ and does not 
affect the seventh place for more than a century. 

15. Since the ecliptic is slowly moving, terms must be added 
to these coefficients, depending on the motion of the ecliptic ; but 
these will be of a higher order still ; and enough has been said to 
show the general character of the coefficients. Returning to the 
general form for them 

£=&!# -h -K a 3 z 

r}=zb x x + b 2 y +b 3 z 
4- c 2 y 4- c. d z 

Suppose (x, y , z) refer to any standard epoch—say 1900*0. 
Then, to find the coordinates for any other epoch, we want the 
nine coefficients, a l9 a 2 . . . c 3 . Now it would not be a difficult 
matter to tabulate these for every year for three hundred years— 
say, from A.D. 1700 to 2000 (all we are likely to want at present)— 
and we should then have the means of bringing up to 1900*0 
accurately any stellar positions referred to another epoch. There 
is no special difficulty about polar stars; the formula? never 
become less simple in any part of the celestial sphere ; and the 
labour of multiplication can be quickly performed with an arith¬ 
mometer, or in many cases by Crelle’s Tables, since the coefficients 
are small. 

16. To reduce from one epoch to another when neither is 
1900*0, we must form nine new coefficients ; but these are readily 
formed when we have those connecting both epochs with 1900*0. 
Thus let (£,17, 0 refer to one epoch, as above, and (X, Y, Z) to 
another, and let 

X=A 1 .'r-{- A 2 y-J- A 3 £, &c. 

SO that £r=A 1 X + B 1 Y-f 0 1 Z, &c. 


according to the usual inversion of direction cosines. 
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Then f =a x x -f a 2 y -f a 3 z 

=: ( <x 1 Ai +^2A2 + <^3A 3 )X + ((x ] B 1 + a 2 B 2 +a 3 B3)Y 

+ (®i C j + a 2^2+« 3 C 3 )Z 

with similar equations for rj and 4, and the law of formation of 
the coefficients is obvious. 

17. Aberration .—The use of direction cosines does not always 
simplify the formulae, however, and we may take aberration as an 
instance where there is no great simplification. The effect of 
aberration is well known to be precisely similar to that of 
parallax, viz. a displacement of the observer in a given direction. 
If the direction be defined by the cosines (l, m, n), and k be the 
coefficient of aberration, the centre of the sphere is virtually 
moved to the point (— kl, — km, — kn). Thus the coordinates of 
a point (; x , y, z) become 

p(x j-kl), p(y + km), p(z + kn), 

p being a constant introduced in order that the sum of the squares 
may remain unity ; i.e. 

p 2 + 2p 2 k(ln+my + nz) -\-p 2 k 2 =1 

How &=2o"*5 and & 2 =o' /, oo2, which is insensible in practice. 
Thus neglecting h 2 , we have 

p—\ —h cos 6 


where cos 0 =lx + my-\-nz, so that 0 is the “Earth’s Way,” 
the angle between the direction of the star and that of the 
Earth’s motion ; and to the first order of k we have for the incre- 
iments of x, y, z due to aberration 


h(l—x cos 0), h{m—y cos 0), h{n—z cos 0). 

Thus the increment of any coordinate x due to precession, 
mutation, and aberration, is 

kl+a l x + a 2 y-^a^z — kx(lx + my+nz), 


^and as regards the calculation of star-corrections, we do not 
appear to gain anything by the use of direction cosines. 

18. It may be remarked incidentally that the effect of aberra¬ 
tion on “ standard coordinates ” 


£=- and 

z z 


is to change them into 


and , 

z-\-hn z+kn 
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the common factory disappearing in the division and leaving only 
the constants (for all stars) hi, km, hn. When we are dealing 
with a small photographic plate for which % is nearly unity, we 
may put z+kn^z^i + hn), and hl/{z-jrhn)=hl, 
and thus £'=£(i— hn) + hl. 

There is thus a displacement of the centre by the quantities 
{Id, h?n) and a change of scale, £ and 77 being both multiplied by 
(1 — hn). This result is derived in another way in Monthly 
Notices, vol. liv. p. 20. 

Summary, 

The present paper points out the advantages of using direction 
cosines or rectangular coordinates in certain astronomical com¬ 
putations, instead of curvilinear coordinates. 

{a). In mapping the surface of the Sun or Moon, observations 
can be easily corrected so as to give us an orthogonal projection 
of the surface on a diametral plane : which is the same as giving 
two of three direction cosines. The third can be easily found, 
and then the transformation to any other axis is made by linear 
formulae. 

(b), In applying precession for long periods to star places, the 
formulae in terms of direction cosines are very simple, and a small 
amount of tabulation would render the accurate reduction of one 
catalogue to another a simple matter. 


The Extra-Equatorial Currents of Jupiter in 1899. By 
Rev. T. E. R. Phillips. 

Despite the southern declination of Jupiter in 189 9, an immense 
amount of detail was visible on the planet’s surface, and numerous 
observations have been received from several observers, which 
have enabled the rotation periods of the surface material in 
various latitudes to be determined with, it is believed, consider¬ 
able accuracy. The writer secured a large number of transits of 
spots during the earlier and intermediate months of the appa¬ 
rition. Mr. Denning, most fortunately, was able to pursue his 
observations very late—in some cases up to the middle of 
September; Mr. A. S. Williams has forwarded a large and 
valuable series of transits; and, in addition to these, much 
assistance has been derived from the figures of Messrs. Gledhill, 
Antoniadi, and J. Comas Sola, of Catala, Spain (the observations 
of the latter being published in Astronomische Nachrichten , No. 
3596, Band 150). Many of the spots have thus been followed 
for considerable periods of time (in some cases for more than 
seven months), and in the majority of cases the observations are 
so numerous and so accordant as to make the question of identi- 
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